Ordered and self— disordered dynamics of holes and defects in the one— dimensional 

complex Ginzburg— Landau equation 
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We study the dynamics of holes and defects in the ID complex Ginzburg-Landau equation in 
ordered and chaotic cases. Ordered hole-defect dynamics occurs when an unstable hole invades a 
plane wave state and periodically nucleates defects from which new holes are born. The results of a 
detailed numerical study of these periodic states are incorporated into a simple analytic description 
of isolated "edge" holes. Extending this description, we obtain a minimal model for general hole- 
defect dynamics. We show that interactions between the holes and a self-disordered background are 
essential for the occurrence of spatiotemporal chaos in hole-defect states. 
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The formation of local structures and the occurrence 
of spatiotemporal chaos are the most striking features 
of pattern forming systems. The complex Ginzburg- 
Landau equation (CGLE) 



A t = A + (1 + i Cl )V 2 A - (1 - ic 3 )\A\ 2 A 



(1) 



provides a particularly rich example of these phenomena. 
The CGLE describes pattern formation near a Hopf bi- 
furcation and has become a paradigmatic model for the 
study of spatiotemporal chaos ■ Defects occur when 
A goes through zero and the complex phase ^:=arg(.A) 
is no longer defined. In two and higher dimensions, such 
defects can only disappear via collisions with other de- 
fects, and act as long-living seeds for local structures like 
spirals || and scroll waves (J] whose instabilities lead to 
various chaotic states For the ID CGLE, however, 

defects occur only at isolated points in space-time (see 
Fig. |l|) and intricate dynamics of defects and local hole 
structures occurs, especially in the so-called intermittent 
and bi-chaotic regimes ||. The holes are characterized 
by a local concentration of phase-gradient q:= d x ip and 
a depression of \A\ (hence the name "hole"), and dynam- 
ically connect the defects (Fig. |l|) . We divide these holes 
into two categories: coherent and incoherent structures. 

Coherent structures - By this we mean uniformly prop- 
agating structures of the form A(x, t) = e~ %u)t A(x — vt) 
[pi . Recently, hole solutions of this form called homoclinic 
holes were obtained || . Asymptotically, homoclinic holes 
connect identical plane waves where A ~ e i(9exx-ut) _ 
With ci, C3 and q ex fixed, unique left moving and unique 
right moving coherent holes are found. Left (right) mov- 
ing holes with q ex = Q (q ex = —Q) are related by the 
left-right q<->—q symmetry of the CGLE. Coherent holes 
have one unstable core mode Q . 

Incoherent structures - In full dynamic states of the 
CGLE, one does not observe the unstable coherent ho- 
moclinic holes, unless one fine-tunes the initial conditions 
(see Fig. Sd). Instead evolving incoherent holes that can 
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FIG. 1. (a) A space-time grey-scale plot of \A\ (dark: 
A « 0), showing the propagation of incoherent holes into a 
plane wave state. The dark dots correspond to defects. Note 
the roughly constant velocities at which the holes propagate. 
Parameter values are ci = 0.6, C3 = 1.4 , with an initial condi- 
tion given by Eq. (@), with 7 = 1, q ex = — 0.03. This non-zero 
q ex breaks the left -right symmetry and results in the differ- 
ing periods of the left and right moving edge holes, (b-c) 
Close-up of I .A I (b) and the complex phase ifi (c). 

grow out to defects occur (Fig. |l| and ||b). 

In this Letter we study the hole — ► defect and defect 
—►holes dynamical processes of the ID CGLE We 
present a minimal model for hole-defect dynamics that 
describes the full "interior" spatiotemporal chaotic states 
of Fig. |l|a, where holes propagate into a self-disordered 
background. Similar "self-replicating" patterns are ob- 
served in many other situations, e.g., reaction-diffusion 
models j]Io| , film-drag , eutectic growth Jl2| , forced 
CGLE ||l3tl and space-time intermittency models |p^| . 

Hole — > defect - Let us consider the short-time evolu- 
tion of an isolated hole propagating into a plane wave 
state. Holes can be seeded from initial conditions like: 



A = exp(i[q ex x + (ir/2) tanh(7a;)]) 



(2) 



The precise form of the initial condition is not impor- 
tant here as long as we have a one-parameter family of 
localized phase-gradient peaks. This is because the left 
moving and right moving coherent holes for fixed c\ , C3 
and q ex are each unique and have one unstable mode only. 
As 7 is varied three possibilities can arise for the time 
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FIG. 2. (a) Schematic representation of the phase space 
of the CGLE around the homoclinic hole solution, showing: 
the ID unstable manifold W u ; the high dimensional neu- 
tral/stable manifold W na that separates decaying from de- 
fect forming states; the manifold T representing the family of 
peaked initial conditions of the form (Q). (b-d) Four snap- 
shots (At = 10) of the g-profile of a right moving hole where 
qex=0, ci=0.6, and C3 = 1.4. The peaked initial condition is 
given by Eq. (^): (b) A hole evolving to a defect (7 = 0.568), 
(c) a decaying hole (7 = 0.5), and (d) a hole evolving close to 
a coherent structure (7 = 0.5545). 

evolution of the initial peak: evolution towards a defect 
(as in Fig. 0a), decay, or evolution arbitrary close to a 
coherent homoclinic hole (see Fig. ||). 

The hole propagation velocities are much larger than 
the typical group velocities in the plane wave states: the 
holes are thus only sensitive to the leading wave. Their 
internal, slow dynamics determines their trailing wave. 
A (nearly) coherent hole will, due to phase conservation, 
have a trailing wave (nearly) identical to the leading wave 
(Fig. H); hence the relevance of the homoclinic holes. 

Defect —> holes - What dynamics occurs after a defect 
has been formed? A study of the spatial defect profiles 
reveals that they consist of a negative and positive phase- 
gradient peak in close proximity (the early stage of the 
formation of these two peaks can be seen in Fig. Ifa; see 
also Fig. 4d of ||). The negative (positive) phase gradi- 
ent peak generates a left (right) moving hole. The life- 
times of these holes depend on their parent defect profile 
(analogous to what we described in Fig. ||) and also on 
ci , C3 and q ex . Hence the defects act as seeds for the 
generation of daughter holes (see also Fig. |l|) . 

Periodic hole-deject states - When an incoherent hole 
invades a plane wave state and generates defects, stable 
periodic hole — > defect — > hole behavior can set in at the 
edges of the resulting pattern jL5| (Fig |l|a) . The asymp- 
totic period r of this process depends on c±, C3, the prop- 
agation direction and the wavenumber q ex of the initial 
condition only; we focus here on right moving holes. The 
period r diverges at a well-defined value of q ex = q co h 
(Fig. |a). This can be understood in the phase space 
picture presented in Fig. |^. Suppose we fix ci and C3. 
The edge defects that are generated periodically yield 
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FIG. 3. (a) Log-linear plot of the period r as a function of 
qex — q C oh- (b) 100/g m as a function of the time At before the 
formation of a defect; (c) q m as a function of q m . 

constant initial conditions for their daughter edge holes, 
similar to fixing 7 in Eq. (^]). The period r will depend 
on the location of the defect profile with respect to the 
stable manifold of the coherent hole. When q ex is varied, 
both this manifold and the defect profile may change, 
and for a certain value of q ex which we call q co h, the de- 
fect generates an initial condition precisely on the stable 
manifold of the coherent hole. The lifetime of the result- 
ing daughter hole then diverges (see Fig. ||d). 

To substantiate this intuitive picture, we have per- 
formed numerics on the dynamics of "edge-holes" invad- 
ing a plane wave state where A ~ e i{q^x-^t) _ v^g have 
performed runs for many different parameters, but will 
only discuss a representative subset here. Our results 
indicate that the r divergence is of the form 



-sln(<7 e 



q co h 



)+TQ 



(3) 



This equation, and in particular the value of s can be 
understood by considering the flow near the saddle point 
shown in Fig ^a. Just after the hole has been formed, 
it first evolves rapidly along the stable manifold. Sec- 
ondly it evolves slowly along the unstable manifold be- 
fore being shot away towards the next defect. For val- 
ues of q ex close to q co h , the holes approach the coherent 
structure fixed point very closely, and r will be domi- 
nated by a regime of exponential growth close to this 
fixed point. Small changes in q ex will have a negligible 
effect on the duration of the first phase (to), but the du- 
ration of the second phase will diverge logarithmically as 
— (1/A) hi(q ex — q co h)- Here A, which depends on c\ and 
C3, denotes the unstable eigenvalue of the coherent struc- 
tures at q ex = q co h ■ In Table 1 we list some numerically 
determined values for q co h, 1/A, and s. We obtained s 
and q co h from a fit of r to Eq. (|3[), whereas A is obtained 
from a shooting algorithm, see Ref. Q. The agreement 
between s and 1/A is quite satisfactory. 

We will now construct a phenomenological model for 
isolated incoherent holes, (i) We will ignore their early 
time attraction to the unstable manifold, and think of 
their location on W v as an internal degree of freedom, 
parameterized by the phase-gradient extremum q m . (ii) 
Clearly the model should have an unstable fixed point for 
values of q m corresponding to coherent holes. We have 
found that, in good approximation, coherent holes have 
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Table 1. Comparison of 1/A with s (see text for details). 

+ gq ex where q n denotes the value of q m for a 
coherent hole in a q ex = state, and g is a negative phe- 
nomenological constant. (Hi) When approaching a de- 
fect, q rn diverges as (At)^ 1 |l6| ]; we have confirmed this 
by accurate numerics (Fig. |3|b) . An appropriate equation 
incorporating these three features is 

qm = A(q m - (q n + gq ex )) + /-t(<7m - {q n + gqex)) 2 , (4) 

where g and /i are phenomenological constants. The first 
term on the RHS of (|]) results from the linearization 
near the coherent fixed point. Nonlinear terms of higher 
than quadratic order on the RHS of Eq. ([|) are ruled out 
by the (At)^ 1 divergence of q m . Our numerical data for 
q m versus q m indeed shows quadratic behavior for large 
enough values of q m (Fig. ||c). For smaller values of q m , 
the curves are quite intricate; this corresponds to the 
rapid early time evolution along the stable manifold not 
included in model ([IJ) . From Eq. <Q) , it is straightforward 
to show that the hole lifetime r (the time taken for q m 
to diverge) displays the required logarithmic divergence 
as q ex is tuned towards a critical value q co h- 

Disordered dynamics - If the patches away from 
the holes/defects were simply plane waves with fixed 
wavenumber, then one would expect, following the argu- 
ments given above, quite regular dynamics. The coupling 
between holes and the background induced by phase con- 
servation becomes the key ingredient to understand dis- 
order in hole-defect dynamics such as shown in Fig [l]a. 
Let us introduce a variable 4> '■— J dxq that measures the 
phasedifference across a certain interval. 

Consider again an edge hole evolving towards a defect. 
While the peak of the g-profile grows, the hole creates 
a dip in its wake (see Fig. ||b) in order to locally con- 
serve 4>. Clearly the trailing edge of this incoherent hole 
is not a perfect plane wave. In the interior of states such 
as shown in Fig. [l]a, unstable holes move back and forth 
through a background of disordered q ex and amplify this 
disorder. Nevertheless, as wc pointed out earlier, the dis- 
ordering dynamics is sufficiently slow such that the holes 
remain approximately homoclinic for much of their lives. 
Although the typical range of values for the disordered 
q ex is small, the hole lifetimes depend on it sensitively. 
Hence the variation in q ex and <p l& sufficient to explain 
the varying lifetimes found in the interior states such as 
that shown in Fig |l|a. Thus the essence of the spatiotem- 
poral chaotic states here lies in the propagation of unsta- 
ble local structures in a self-disordered background. 



Minimal model - To illustrate our picture of self- 
disordered dynamics, we will now combine the various 
hole-defect properties with the left -right symmetry and 
local phase conservation of the CGLE to form a mini- 
mal model of hole-defect dynamics. From our previous 
analysis, we see that the following hole-defect properties 
should be incorporated: (i) Incoherent holes propagate 
either left or right with essentially constant velocity (see 
Fig. |l|a). (ii) For fixed ci,C3, their lifetime depends on 
the profile of their parent defect, the direction of prop- 
agation, and on the wavenumber of the state into which 
they propagate. (Hi) Eq. ([|) captures essentially all as- 
pects of the evolution of their internal degree of freedom. 
When q m diverges, a defect occurs. 

In our model we will assume that all the defects have 
the same profile and so act as unique initial conditions 
for their daughter incoherent holes. While in principle a 
defect profile could depend on the entire history of the 
hole which preceded it, for simplicity we have chosen to 
neglect this. We have observed that for some regions of 
the ci, C3 parameter space, the defect profiles from the in- 
terior spatiotemporal chaotic patterns show a surprising 
lack of scatter |J . Therefore we believe that treating the 
defect profiles as constant, and only including the effect 
of the background in the hole dynamics incorporates the 
essence of the coupling to a disordered background. 

We discretize both space and time by coarse-graining, 
and take a "staggered" type of update rule which is com- 
pletely specified by the dynamics of a 2 x 2 cell (see 
Fig. |]a). We put a single variable <f>i on each site, 
corresponding to the phase difference across a cell di- 
vided by 27r. Local phase conservation is implemented 
by + 4>' r — 4>i +</v, where the primed (unprimed) vari- 
ables refer to values after (before) an update. 

Holes are represented by active sites where |</>| > 4>t] 
here </> plays the role of the internal degree of freedom. 
Inactive sites are those with \<p\ < (j) tl and they repre- 
sent the background. The value of the cutoff <fi t is not 
very important as long as it is much smaller than typical 
values of for coherent holes. Here (f>t is fixed at 0.15. 
Without loss of generality we force holes with positive 
(negative) <j) to propagate only from 4n (<f> r ) to 4>' r (4>'i). 

Depending on the two incoming states, we have the 
following three possibilities: 

One site active: Without loss of generality we assume 
that we have a right moving hole. We implement evolu- 
tion similar to Eq. (Q) , but neglect the quadratic term of 
Eq. (JiJ); even though q m diverges, the local phasediffer- 
ence 4> m does not diverge near a defect. Hence the finite 
time divergence of the local phase gradient q that signals 
a defect can be replaced by a cutoff <fid for <f>. Therefore, 
when (f>i < 4>d, the internal hole coordinate <p is taken to 
evolve via <j)' r = 4>i +A(0; — <fi n — g4> r )- Here A sets the 
time scales and can be taken small (fixed at 0.1). This 
evolution equation, combined with the local phase con- 
servation, means that an incoherent hole propagating 
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FIG. 4. (a) Grid model geometry showing the sites (dots) 
and hole propagation direction (arrows). The update rule is 
defined within a 2 x 2 cell, (b-d) Dynamical states in the grid 
model, for 4>„ = 0.6 and 4> a d = 0.75. Initial condition: center 
site has = 0.7, everywhere else <f> = (hence the symmetric 
patterns), (b) g — and 0<j = l. (c) Disordered dynamics for 
nonzero coupling (g— —3, <f)d = 1). (d) Zigzag structures occur 
for 3 = -3, (j> d = 0.98. 

into a perfect laminar state will leave a disordered state 
in its wake. When <pi > cf>d, & defect occurs and two new 
holes are generated: 4>' r = </w, and §\ — <f>d — 1 — (fad- 
The factor —1 reflects the change in winding number at 
a defect. 

Both sites inactive: Away from the holes/defects, the 
relevant dynamics is phase diffusion. This is implemented 
via: (j)' r = D(f)i+{l-D)(f) r . The value of D is fixed at 0.05 
and is not very important. 

Both sites active: This corresponds to the collision of 
two oppositely moving holes. Typically this leads to the 
annihilation of both holes (see Fig. |l|a), which we imple- 
ment here via phase conservation: (f)' r =</>[ = {4>i+4>r) /2. 

The coupling of the holes to their background, g, 
should be taken negative (although its precise value is 
unimportant). For g — the lifetime r becomes a con- 
stant, independent of the <j> of the state into which the 
holes propagate, and moreover, the dynamical states are 
regular Sierpinsky gaskets (Fig. [|b). Nevertheless, start- 
ing from a 4> — state, the local phase conservation of the 
hole dynamics leads to a background state with a disor- 
dered <[> profile. For g < the coupling to this background 
leads to disorder as shown in Fig. [|c,d. This illustrates 
the crucial importance of the coupling between the holes 
and the self-disordered background. 

The essential parameters determining the qualitative 
nature of the overall state are <p n , 4>d and 4> a d- These pa- 
rameters determine the amount of phase winding in the 
core of the q ex = coherent holes (<j> n ) and in the new 
holes generated by defects (4> a d,4>d — 1 — 4>ad)- When 
varying the CGLE coefficients ci,C3, these parameters 
change too; for example, (f) n typically decreases when c\ 
or C3 are increased. As a result, for large values of c\ and 



C3, and <j)' r are typically larger than <f> n so that most 
"daughter holes" will grow out to form defects and hole- 
defect chaos spreads (Fig. 4c, d). For sufficiently small 
values of c\ and C3, on the other hand, <fi„ is large and 
both daughter holes will decay. For intermediate values 
of c\ and C3 it may occur that is significantly larger 
than cf)' r , leading to zigzag states [6] (Fig. 4d). 

In conclusion, we have studied in detail the dynamics 
of local structures in the ID CGLE. We have obtained a 
quantitative understanding of the edge holes, unraveled 
the interplay between defects and holes, and put forward 
a simple model for some of the spatiotemporal chaotic 
states occurring in the CGLE. 
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